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Abstract 

There has been observed widespread destruction of natural ecosystems around the world due to population growth, 
land use change and clear cutting which have affected soil properties. Different management strategies have been so 
far implemented to reduce this crisis in various regions of the world, such as e.g. short-term and long-term conserva- 
tion management in the Zagros region. However, any management approach should be evaluated with appropriate 
measures to determine how managed areas respond. The main objective of the present study was to evaluate the 
potential of earthworms as an indicator for different forest management strategies and human disturbances in Zagros 
oak (Quercus persica Jaub. and Spach) forest. The sites selected included undisturbed one as the control (Un), the 
sites under five-year conservation management (FCM) and twenty-year conservation management (TCM) as well as 
the disturbed site (D). The results of principal component analysis (PCA) showed that different regions separated into 
the components: PCI and PC2. Un and TCM sites gathered together and represented higher values of the factors such 
as pH, K available , OC, clay content, P ava ;i a bi e , CEC, overstory tree canopy, N tot , biomass and abundance of earthworms. 
The positive direction of the first axis reflected a gradient of EC, BD and P tot . According to the logistic model, NH 4 -N 
and EC played the most important role in earthworm presence and absence in Zagros forest ecosystem. Earthworm 
abundance and biomass could be a good indicator to evaluate different forest management strategies in the study area. 
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Introduction 

There is a strong relationship between human economic 
activities and natural resources, especially in forests 
and pastures in the Zagros province. Man’s activities 
have caused different disturbances including soil ero- 
sion and clear-cut forest areas in addition to decline of 
soil fertility as well as species diversity and richness 



(Salehi et al. 2011). In the region,, different management 
strategies have been implemented, such as e.g. short- 
term and long-term conservation management, which 
affected its ecosystems in different ways (Zandebasiri 
et al. 2010). Soil quality is one of the prominent factors 
in the healthiness of different ecosystems. The knowl- 
edge on the processes ongoing in the ecosystem as well 
as that on soil as a vital part in the ecosystem is essential 
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for finding good solutions, and therefore the properties 
of soil which are appropriate for evaluation of ecosystem 
functions should be better known and defined (Karlen 
et al. 2003). To study soil quality scientists have offered 
different indicators such as physical (Hseu et al. 1999) 
and chemical and biological soil properties (Doran and 
Parkin 1994). For monitoring soil quality indices, it is 
better to select soil properties that change with chang- 
ing management approaches, and also show good cor- 
relation with the processes occurring in soil (Herrick 
2000). Studies on soil quality with the aim to identify 
the effects of different management types on forests de- 
struction are imperative (Khaledian et al. 2011). Hence, 
if their results point to the effects of management on soil 
quality, they can constitute suitable support in achiev- 
ing better understanding of sustainable management 
(Herrick 2000; Bhardwaj et al. 2011). 

Soil organisms play a vital role both in sustainable 
management and maintaining soil quality as they in- 
fluence various properties of soil on the one hand, and 
on the other- are influenced by soil nature (Brady and 
Weil 2008). One of the most important issues for a pro- 
grammer is the question how changes in soil biological 
communities influence soil quality (Chan and Barchia 
2007). Most studies address the impacts of management 
strategies and destructions due to human activities on 
forests ecosystems through monitoring of forest plant 
diversity and richness (Kouame et al. 2004), while there 
remain scarce knowledge about soil organisms such as 
earthworms. There is a wide destruction of forest eco- 
systems around the world owing to population growth, 
land use change and forest clear cutting which have 
strongly affected soil fertility (Dinesh et al. 2003) and 
soil biology (Islam and Weil 2000). Investigations show 
that some soil biological characteristics, such as changes 
in soil biological communities are measurable, in other 
words a study on soil microorganisms and effective pro- 
cesses on them can be a key indicator of soil quality 
in the future (Acost-Martines et al. 2007). Doran and 
Zeiss (2000) state that the evaluation of soil organisms 
as a good indicator of soil quality can be more inves- 
tigable, and also less time-consuming and less costly. 

Among soil organisms, earthworms are the most 
commonly occurring invertebrates (Tondoh et al. 
2007). Earthworms enhance soil fertility and produc- 
tivity through increasing availability of mineral nutri- 
ents to plants (Brady and Weil 2008) and improving 



the dynamics of soil organic matter (Reich et al. 2005). 
Abundance and activity of earthworms are influenced 
by soil conditions. Mboukou- Kimbasta et al. (2007) de- 
clared that soil nitrogen had a positive effect on earth- 
worm abundance and biomass. Earthworms mineralize 
phosphor, nitrogen, and organic matter in manure ferti- 
lizers and turn them into nutrients available for plants 
(Buchanan et al. 1988). These organisms directly in- 
crease litterfall decomposition rate by enhancement of 
microbial activities (Tsukamoto and Sabang 2005) and 
advance soil physical properties by improvement of soil 
structure as well as water and air infiltration capabil- 
ity which are indispensable for plant growth (Sautter 
et al. 2006). The effect of land use change on density 
and abundance of earthworms in Hariyali area in India 
was examined by Bhadauria et al. (2012). The results 
indicated the lowest abundance of earthworms in natu- 
ral forest which was severely destructed. The authors 
showed that conservation practices including restora- 
tion of destroyed forest improved earthworm commu- 
nities and organic matter for the duration of 20 years. 
Tondoh et al. (2011) investigated the effects of forest de- 
struction on earthworm population in the western part 
of Ivory Coast. The results showed that forest distur- 
bance decreased soil organic carbon and soil pH. At the 
same time, multi-variable regression analysis between 
earthworm communities and environmental variables 
showed that soil organic matter and pH were potential 
indicators to reflect earthworm abundance. 

There are studies about changes in abundance, pres- 
ence or absence and biomass of earthworms observed 
under destructive activities or diverse management types 
in different regions (Hendrix and Bonlen 2002; Aubert et 
al. 2003; Haynes et al. 2003; Rombke et al. 2005 ; Ammer 
et al. 2006). Yet, previous studies have seldom been fo- 
cused on both soil physiochemical properties and earth- 
worms under destruction and short- or long - term con- 
servation management conditions. Regarding this topic, 
no comprehensive study in Zagros oak forests has been 
to this point carried out. The present study has been per- 
formed for the first time in Zagros oak forests located in 
the western part of Iran. The results of the study can offer 
managers and researchers novel management approaches 
to examine the potential of earthworms as a biological 
indicator for evaluation of predominating conditions in 
forest ecosystems. The main objective of the study was 
to evaluate the potential of earthworms as an indicator 
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for different forest management strategies and human 
disturbances in the ecosystem. 

Material and methods 

Study area 

The study area is located in Zagros forests in the western 
part of Iran (fig. 1). In 2012, we chose four sites within 
100 ha of forest area, The sites selected were divided 
into undisturbed one as the control (Un), the sites under 
five-year conservation management (FCM) and twenty- 
year conservation management (TCM) and the disturbed 
site (D). Prior to the implementation of management prac- 
tices, the conditions of FCM and TCM sites were simi- 
lar to those at D site. According to an aerial photograph 
taken in 1965 as well as communication with experts and 
local community, all the sites were then covered by dense 
forest. All the sites chosen are alike in terms of physi- 
ography (slope and altitude) and tree species structure. 



The dominant tree species is oak ( Quercus persica Jaub. 
and Spach) which covers more than 90% of the study 
area. Average annual precipitation and temperature are 
590 mm and 17°C, respectively. 

Soil sampling 

Soil samples were taken from 0-30 cm soil layer in May 
2012, using the systematic-random sampling method 
(100 m x 200 m). Further analyses of soil physical and 
chemical properties were carried out on composite sam- 
ples consisting of three sub-samples taken from select- 
ed soil layer (Maranon et al. 1999). The soil samples 
were air dried and sieved (2 mm) to remove roots and 
coarse fragments. 

Analytical methods 

Soil texture was determined by the Bouyoucos hydrom- 
eter method (Bouyoucos 1927). Soil bulk density (BD) 
was determined by the undisturbed soil core method 
(Blake and Hartage 1986). Saturation percentage was 





Folia Forestalia Polonica, series A, 2014, Vo I. 56 (3), 121-134 



124 



Mehdi Heydari, Hassan Poorbabaei, Masoud Bazgir, Ali Salehi, Javad Eshaghirad 



determined gravimetrically (Famiglietti et al. 1998). 
Soil pH and electric conductivity (EC) were measured 
using a glass -electrode pH meter and EC meter, respec- 
tively (Kalra and Maynard 1991). Organic carbon (OC) 
was determined based on the Walkely and Black rapid 
titration method (Black 1979). Total nitrogen (N) was 
measured by the Kjeldah method (Bremmer and Mul- 
vaney 1982). Flame photometry (Bastida et al. 2007) 
was used for measurement of total potassium. Soil 
cation exchange capacity (CEC) was determined based 
on NH40Ac extraction (Sumner and Miller 1996). Ex- 
changeable cations (Al 3 h , Fe 3+ , Mn 2+ , Zn 2+ , Ca 2+ , Mg 2+ , 
K + , Na + ) were extracted using a 0.1 M BaCl 2 (Hender- 
shot and Duquette 1986). The extracts obtained were 
centrifuged, filtered (Machery-Nagel 615 1/4 filters 
rinsed with BaCF) and analysed by ICP-AESS (Varian, 
Australia). Available phosphorus was assessed by the 
Bray and Kurtz method (1945), soil N0 3 -N and NH 4 -N 
concentrations (extracted using 2M KC1 solution) were 
measured with a spectrophotometer (Maynard and Kal- 
ra 1993). Total P was measured calorimetrically using 
the ascorbic acid method, followed by digestion in nitric 
and perchloric acid (Sommers and Nelson 1997). 

Earthworm sampling 

Earthworms were collected manually, around the center 
of each sample plot- in the lateral soil area of 50 x 50 cm 
drilled down to the depth of 20 cm. Earthworms were 
counted and then dried in the oven at 60°C for 24 h. 
Next, dried earthworms were weighed to determine bio- 
mass (Welke, Parkinson 2003). Earthworm abundance 
was calculated as the number of earthworms per surface 
area (m 2 ). 

Statistical analyses 

The normality of data distribution and homogeneity of 
variances were verified by Kolmogrov- Smirnov and 
Levene’s tests, respectively Environmental data were 
analyzed using one-way ANOVA and in line with ho- 
mogeneity of variance, significant differences between 
means were determined by Duncan’s test. Consistent 
with data normality we used Pearson’s correlation to 
analyze the relation between earthworm densities (No/ 
m 2 ) and biomass and environmental variables. In order 
to determine the most effective factors with influence 
on plot distribution, multivariate procedure such as 
principal component analysis (PCA) was used (Arekhi 



et al. 2010). PCA was used in order to trim down many 
variables to a smaller group of newly derived variables 
(components) that adequately summarized initial infor- 
mation. These variables were independent of each other, 
and for that reason could be used as predictors in PCA. 
PCA is the ordination technique that constructs a theo- 
retical variable that minimizes the total residual sum of 
squares after fitting straight lines to species data. Data 
standardization was necessary to use PCA because the 
variables were measured in different units. Therefore, 
data were centered and standardized by the standard de- 
viation (Jalilvand and Kooch 2012). Logistic regression 
was used to investigate the relationship between earth- 
worm presence and absence and environmental factors. 
Also, the power and efficiency of the logistic regression 
model was evaluated based on the omnibus test. Two 
statistics coefficient of determination’s Pseudo (Cox and 
Snell R Square and Nagelkerke R Square) were used 
to examine the explanatory power of independent vari- 
ables in expression of a dependent variable (Habibpoor 
and Safri, 2012). The coefficient of determination’s 
Pseudo is between 0 and 1. If the value of this statistic 
is closer to 1, it means that independent variables have 
a greater role in explanation of variance in a dependent 
variable. In addition, the Hosmer-Lemeshow test was 
used to determine fitness of the logistic model. To run 
the model, the following equation was used: 



log it(Y) = natural log (odds) = In 



P 



1 -P 



— a + /), X 1 + $ 2 X 2 + Pk^k 



where: 

P - the probability of an event with the presence of an 
independent variable X, 

a - the constant or parameters related to the coordinate 
axes, 

ft - the regression coefficient and X represents inde- 
pendent or predictive variables. 



Simple linear regression was used to predict the val- 
ue of a dependent variable (y) (here the most important 
soil attributes in the presence or absence of earthworms, 
based on logistic regression) derived from the value of an 
independent variable (x) (here biomass or abundance of 
earthworms) based on the following equation: 

y = a + bx 
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where: 

y - a dependent variable, 
a - the intercept, 
b - a slope, 

x - an independent variable. 

The statistical analyses were carried out using PC- 
ORD, V. 4.17 package (McCune and Mefford 1999) and 
SPSS, versionl6. 

Results 

Principal component analysis (PCA) 

To distinct the most effective variables on the distribu- 
tion of sample plots, PCA was performed for 17 varia- 
bles in 200 plots. The results of PCA ordination are pre- 
sented in tab. 1 and fig. 2. Broken-stick eigenvalues for 
the data set indicated that the first two principal com- 
ponents (PCI and PC2) clearly captured more variance. 
The first two principal components together accounted 
for 63% of the total variance in the data set. Therefore, 
54.38% and 8.64% variances were accounted for the 
first and second principal components, respectively 
(tab. 2). The first principal component was by far the 
most important for representing variation in this study. 



Table 1. Pearson’s correlation coefficient of environmental 
variables and PCI and PC2 



Variables 


AXIS 1 


AXIS 2 


1 


2 


3 


pH 


-0.6** 


0.236 ns 


EC/dS nr 1 


0.868" 


0.13 ns 


Pava/mg kg 1 


-0.898" 


0.147 ns 


Kava/mg kg- 1 


-0.696” 


0.07 ns 


%/Clay 


-0.887” 


0.105 ns 


%/Silt 


0.168 ns 


-0.711” 


%/Sand 


0.278 ns 


0.155 ns 


BD/g cm 3 


0.773** 


0.444* 


%/OC 


-0.903” 


0.167 ns 


%/N tot 


-0.922** 


0.171 ns 


P to/mg kg- 1 


0.852” 


0.246 ns 


K, ot /mg kg- 1 


-0.304* 


-0.462* 


N03-N/mg kg- 1 


-0.07 ns 


-0.832” 


NH4-N/mg kg- 1 


-0.931” 


0.130 ns 



1 


2 


3 


CEC/cmol (+) kg -1 


-0.818** 


0.139 ns 


Overstory tree canopy 1% 


1 

© 


-0.09 ns 


Earthworms/No nr 


-0.551* 


-0.23 ns 


Earth worm biomass gr/m 2 


-0.652** 


-0.229 ns 



* Significant (P < 0.05),** Significant (P < 0.01), ns: No 
significant ( n = 50). 




Axis 1 



O Undisturbed or control site (Un) 

V Disturbed site (D) 

□ 20 years conservation management site (TCM) 

O 5 years conservation management site (FCM) 

Figure 2. PCA ordination diagram for sample plots of 
different sites 

The results show that negative direction of the first 
axis represents the factors such as pH, K available oc clay 
Pavaiiabie? CEC, overstory tree canopy, N tot , as well as 
earthworm abundance and biomass. In other words, 
these factors showed a strong correlation with PCI 
negative direction. Un and TCM sites were strongly 
correlated with the negative direction of the first axis 
and were gathered together. This means that environ- 
mental factors in Un and TCM sites were similar. The 
positive direction of the first axis reflected a gradient 
of EC, BD and P tot . The sample plots in FCM and D 
sites gathered together around this axis (fig. 2). As 
mentioned above, Un and TCM sites indicated higher 
earthworm density. Earthworms showed a positive 
response to increased environmental factors in these 
sites. 
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Table 2. Eigenvalue, variance percentage, cumulative 
percentage of variance and broken-stick eigenvalue from 1 
to 10 axes PCA axes 



Axes 

No 


Eigen- 

value 


Variance 

percentage 


Cumulative 
percentage 
of variance 


Broken-stick 

Eigenvalue 


1 


11.965 


54.385 


54.385 


3.691 


2 


1.902 


8.645 


63.030 


2.691 


3 


1.563 


7.107 


70.137 


2.191 


4 


1.078 


4.900 


75.037 


1.857 


5 


0.964 


4.383 


79.420 


1.607 


6 


0.755 


3.430 


82.850 


1.407 


7 


0.675 


3.070 


85.919 


1.241 


8 


0.467 


2.122 


88.041 


1.098 


9 


0.430 


1.956 


89.997 


0.973 


10 


0.358 


1.629 


91.626 


0.862 



Logistic regression model 

The results showed that the presence and absence of 
earthworms can be explained with high probability 
(73%) by independent variables including soil physi- 
cal and chemical properties and overstory tree canopy. 
Between the variables entered into regression analysis, 
OC, N tot , NH 4 -N, N0 3 -N, CEC, BD, EC and overstory 
tree canopy were able to predict variation of a depend- 
ent variable, i.e. the presence -absence of earthworms. 
The evaluation of the power and efficiency of the 
logistic regression model based on the omnibus test 
showed that the model was acceptable and significant 
(P < 0.01). Two statistics coefficient of determination’s 
Pseudo, Cox and Snell R 2 and Nagelkerke R 2 were as 
high as 0.73 and 0.78, respectively, and this showed 
that independent variables had important role in ex- 
plaining variance of earthworm presence and absence 
of. In fact, these variables explained 73.22 and 78.6% 
percentage the variations of the presence and absence 
of earthworms. 

Hosmer-Lemeshow test for goodness of fit 

Based on the Hosmer-Lemeshow test, fitting predic- 
tion variations of a dependent variable was significant. 
In other words, the research model had good fitness 
and independent variables were able to predict high 
proportion of a dependent variable (earthworms). 



The results showed that the variables entered into 
the model: NH 4 -N, N0 3 -N, pH, CEC, OC and K avail _ 
ablc had a positive effect on the probability of the pres- 
ence of earthworms. Exp (B) is the exponentiated B 
coefficient. When Exp (B) > 1, it means that with in- 
creasing values of an independent variable, the pres- 
ence of earthworms is more probable. BD and EC had 
a negative effect on the probability of the presence of 
earthworms. In other words, with decreasing BD and 
EC values, the possibility of earthworm presence in- 
creased. B is the logistic regression coefficient that is 
significant (+ or -) and shows the effects of each inde- 
pendent variable on a dependent variable. The Wald 
statistic provides a significance test for each predictor. 
If the statistic was significant for a variable, this vari- 
able was entered into the model (tab. 3). 

Table 3. Logistic model for earthworm presence and 
absence, exponentiated coefficients or Exp (B), significance 
level (sig), Wald chi-square test (Wald) and logistic 
regression coefficient (B) 



Variables 


Exp(B) 


Sig 


Wald 


B 


EC/dS nr 1 


0.37 


0.03 


3.11 


1.7- 


BD/g cm' 3 


0.23 


0.01 


5.9 


0.34- 


pH 


1.13 


0.04 


0.8 


0.14 


%/OC 


1.14 


0.03 


4/2 


1.15 


P ava /mg kg' 1 


1.17 


0.02 


4.2 


1.12 


K ava /mg kg' 1 


1.23 


0.03 


2 


3.5 


N03-N/mg kg' 1 


1.33 


0.40 


6 


3.2 


NH4-N/mg kg' 1 


1.45 


0.01 


6.1 


6.18 


CEC/cmol (+) kg' 1 


1.17 


0.03 


4.7 


2.3 



The logistic model for earthworm presence-absence in this 
study was explained by the following equation: 



P (earthworm presence and absence) = 

= 8.5 + 6.18 NH 4 -N + 3.2 N0 3 -N + 3.5 K ava +1.12 P ava + 

+ 1.15 OC+ 0.14 pH + 2.3 CEC - 1.7 EC - 0.34 BD 

Based on the above model, increased NH 4 -N, 
N0 3 -N, K available , P available , OC, pH and CEC had a posi- 
tive effect on the presence of earthworms whereas in- 
creased EC and BD had a negative effect on the pres- 
ence and absence of earthworms. 
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Correlation between the abundance and biomass of earthworms 
and environmental factors 

According to Pearson’s correlation earthworm abundance positively corre- 
lated with pH, clay contents, N tot , NH 4 -N, while EC and BD showed a nega- 
tive correlation with earthworm abundance.. Furthermore, earthworm bi- 
omass showed a positive correlation with pH, OC, P avai i a bi e > K available , clay 
contents, N tot , NH+4 and CEC, while it correlated negatively with EC and 
BD (tab. 4). These results corresponded with the logistic model. 

One-Way ANOVA showed that all the factors examined in this study 
significantly differed among the study areas. Duncan’s test showed that pH, 
NH 4 -N, N tot , CEC, OC, clay contents,, canopy percentage and biomass of 
earthworms were the highest in Un and TCM sites and the lowest - in D and 
FCM sites (tab. 5 and fig. 3). 





Figure 3. Abundance and biomass of earthworms in four sites with different 
conditions 

Un - control; FCM - five-year conservation management; TCM - twenty-year 
conservation management; D - disturbed site. 

The values represent mean ± standard error (n = 50), different lowercase letters 
indicate significant differences (P < 0.01). 

According to the logistic model, NH 4 -N and EC had the most important 
role in the presence and absence of earthworms. Therefore, linear regression 
between abundance and biomass of earthworms and these two important 
factors had shown in fig. 4. There was observed a strong positive correlation 
between earthworm abundance and biomass and NH 4 -N. This means the 
two variables (earthworm abundance and biomass and NH 4 -N) varied in 
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Table 5. Soil characteristics, earthworm abundance and biomass and overstory tree canopy in four sites with different 
conditions 



Sites 


Variables 


Un 


TCM 


FCM 


D 


Sig 


pH 


7.46 (±0.01) a 


7.32 (±0.07) b 


7.21 (± 0.01) c 


7.25 (±0.01) c 


0.02* 


EC/dS nr 1 


0.34 (± 0.01) c 


0. 38 (±0.02) c 


0.67 (± 0.01) b 


0.88 (±0.01) a 


0.04* 


%/OC 


4 (±0.14) a 


2.9 (±0.03) b 


1.33 (±0.02) c 


1.23 (±0.06) c 


0.001** 


Pava/mg kg 4 


22.76 (±0.5) b 


26.71 (±0.5) a 


8.6 (±0.39) c 


8.23 (±0.35) c 


0.001** 


K ava /mg kg" 1 


621.26 (±9) a 


488.16 (±8.8) b 


425.77 (±1.97) c 


421.7 (±5) c 


0.000** 


%/Clay 


36.48 (±0.57) a 


38.7 (±0.4) a 


21.95 (±0.67) b 


20.22 (±0.8) b 


0.002” 


%/Silt 


34.7 (±0.52) b 


35.1 (±0.54) b 


41.54 (±0.8) a 


34.48 (±0.87) b 


0.02* 


%/Sand 


28.74 (±0.73) c 


26.14 (±0.4) c 


34.71 (±1) b 


45.28 (±0.8) a 


0.000** 


BD/g cm" 3 


1.18 (±0.01) d 


1.26 (±0.008) c 


1.31 (±0.009) b 


1.59 (±0.009) a 


0.01* 


%/Ntot 


0.59 (±0.01) a 


0.307 (±0.006) b 


0.12 (±0.002) c 


0.14 (±0.007) c 


0.001” 


P to/mg kg- 1 


549.95 (±5.7) d 


610.38 (±7.3) c 


697.38 (±9.7) b 


880.47 (±6.3) a 


0.000” 


Kto/mg kg- 1 


3362.23 (±31.29) b 


3221.93 (±39) c 


3602.93 (±18) a 


3375.75 (±36) a 


0.001* 


N03-N/mg kg' 1 


22.72 (±0.36) b 


29.86 (±0.43) a 


20.91 (±0.39) c 


15.48 (±0.36) d 


0.002” 


NH4-N/mg kg' 1 


0.64 (±0.01 ) a 


0.66 (±0.008) a 


0.18 0.018 b 


0.15 (±0.003) b 


0.01* 


CEC/cmol (+) kg' 1 


2.24 (±0.03) a 


1.88 (±0.03) b 


1.43 (± 0.01) c 


1.45(±0.01) c 


0.002** 


Earthworms/No nr 


6.3 (±0.84) a 


6.1 (±0.8) a 


4 (±0.6) b 


3.84 (±0.5) b 


0.02* 


Earthworms 
biomass g nr 2 


2.54 (±0.35) a 


2.48 (± 0.35) a 


1.42 (±0.22) b 


0.99 (±0.14) b 


0.000” 


Overstory tree 
canopy/% 


40 (±3.32) a 


35 (± 2.97) a 


25 (±2) b 


15 (±1.5) c 


0.000” 



Values are mean ± standard error (n = 50); different lowercase letters represent difference significant (P < 0.05); * significant 
(a = 5%); ** significant (a = 1%) according to Duncan test. 



the same direction (R 2 = 0.641 and 0.608 respectively) 
(fig. 4A, 4B). At the same time, there was found a strong 
negative correlation between earthworm abundance 
and biomass and EC. A negative correlation means that 
two variables (here, EC and earthworm abundance and 
biomass) vary in opposite directions (fig. 4C, 4D). In 
other words, we can say that EC and NH 4 -N values can 
be predictable based on earthworm abundance and bio- 
mass. 

Discussion 

Studies on soil quality and suitable indices expand with 
broadening human knowledge about soil as an essential 
agent in biosphere, which is essential not only for food 



production, but also for preservation of natural environ- 
ment, at a local, regional, and global level (Glantz 1995). 
Taking into consideration a key role of earthworms in 
organic matter cycling and improving the soil structure 
in numerous ecosystems, the recognition of elements 
which control earthworm abundance is definitely nec- 
essary for soil sustainable management (Fonte et al. 
2009). 

Physical and chemical properties of soil 
and canopy percentage in relation 
to management and destruction 

The results of this study showed that environmental 
factors in the control site (Un) and those under twenty- 
year conservation management site (TCM) are consid- 
erably similar. In Un and TCM sites, the values of soil 
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Figure 4. Linear regression between earthworm abundance 
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biomass and EC and NH 4 -N in the study sites 



pH, CEC, K. lvailable and P available , clay contents, OC, total 
and mineral nitrogen, foliar litterfall percentage and 
canopy percentage were higher than in five-year conser- 
vation management site (FCM) as well as in not man- 
aged or highly disturbed (D) sites. On the other hand, 
in FCM and D sites, soil EC, BD and P tot were higher 
than in Un and TCM sites. Consistent with the results, 
the destructive factors such as overgrazing, understory 
agriculture and so on in FCM and D sites had more ef- 
fects on soil physical and chemical properties. Besides, 
edaphic conditions in FCM site, even though it is desig- 
nated for this conservation practice, are not classified as 
suitable. Long-term conservation management positive- 
ly influenced physical and chemical properties of soil, 
so that environmental conditions in TCM were similar 
to those in Un site. Studies have proved that destruction 
affected soil properties, which resulted in forest frag- 



mentation, and finally changed vegetation and soil nu- 
trient contents in forest (Kelemen et al. 2012). The can- 
opy of trees and shrubs in FCM and D sites was almost 
half of that in Un and TCM sites. This can be one of the 
main factors affecting soil quality in the areas studied. 
One study showed that the canopy of blue oak ( Quercus 
douglasii Hook, and Arn.) significantly increased soil 
fertility (Dahlgren et al. 2003). In contrast, clear cutting 
decreased soil nutrient and organic matter content and 
also had a negative effect on soil nutrient cycling (Dahl- 
gren et al. 2003). Soil organic matter is the main source 
of nutrients and soil fertility, and plays an important 
role in advancing soil physico-chemical and biological 
properties (Brady and Weil 2008). Studies have showed 
positive effects of conservation management on soil 
properties (Shaoshan et al. 2008). Shaoshan et al (2008) 
studied the effects of deforestation on soil quality in the 
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Ziwuling region located in northwestern Shaanxi prov- 
ince (China) and showed that with increasing destruc- 
tion, there drastically changed soil physico-chemical 
properties, and particularly - organic matter contents. 
In this study, high bulk density observed in D site can 
be related to low organic matter content (Ahmad et al. 

2007) and soil compaction (Dauda and Samari 2002). 
In Un and TCM sites, high clay percentage in soil es- 
sentially influenced an increase in nutrient contents and 
CEC. High CEC values found in Un and TCM sites, 
were not unexpected, because both clay and organic 
matter contents (Rashidi and Seilsepour 2008), which 
both have been in these regions higher than in two sites 
of (D) and (FCM). Generally, our results showedthat the 
two conservation management types (FCM and TCM) 
and destruction obviously influenced soil physical and 
chemical properties in Zagros forest ecosystem. 

Abundance and biomass of earthworms 
in relation to soil properties 

It was found the positive effect of conservation man- 
agement on biomass and presence of earthworms. 
Earthworm abundance and biomass in Un and TCM 
sites were almost two times higher than in two de- 
structed sites (D and FCM). The study on Irish pastures 
proved that an increase of earthworm abundance was 
connected with a longer period of conservation man- 
agement, which resulted from augmenting earthworm 
community by adding plant debris into soil (Curry et al. 

2008) . Our results signify the negative effect of destruc- 
tion related to soil properties on abundance, biomass 
and presence of earthworms., Other study indicated 
that alteration and destruction of forests impacted soil 
physico-chemical properties, followed by a decrease 
of earthworm density and abundance (Bhadauria et al. 
2012). The results showed a strong positive relationship 
between earthworm abundance and biomass and N tot , 
NH4-N, OC, Pavaiiabie, K available , clay contents and CEC. 
On the other hand, there was a negative correlation ob- 
served between earthworm abundance and biomass and 
EC and BD. In fact, the suitable condition of soil in Un 
and TCM increased abundance and biomass of earth- 
worms. Boyer et al. (1991) found a positive correlation 
between earthworm abundance and soil fertility. Kooch 
et al. (2010) studied the effect of canopy opening in 
beech forest on soil chemical and biological properties 
and found higher earthworm abundance and biomass in 



the areas with closed canopy, which indicated high val- 
ues of soil C/N ratio and CEC. Mboukou-Kimbatsa et 
al. (2007) declared that soil N positively affected abun- 
dance and biomass of earthworms. On the other hand, 
canopy opening in D site, above and beyond its negative 
effects on different soil properties and consequential- 
ly - on earthworms population, could directly reduce 
earthworm abundance by bringing about an increase of 
soil temperature (Nachtergale et al. 2002). According 
to the results of this study, there is a positive correla- 
tion between earthworm abundance and biomass with 
clay percentage in soil and also according to PCA high 
clay percentage, abundance and biomass of earthworms 
in two sites of (Un) and (TCM) have been approved. 
In a study on earthworms in southern Australian fields, 
there was asserted a negative correlation between earth- 
worm population and sand percentage in soil (Bucker- 
field et al. 1997). In addition, Fonte et al. (2009), when 
investigating earthworm population in relation to soil 
organic matter dynamics and management in California 
cropping systems showed a strong positive correlation 
between earthworm abundance and soil clay percentage 
On the other hand, no correlation between earthworm 
abundance and silt and sand, was shown, which corre- 
sponds with our results. 

Presence and absence of earthworms in relation 
to soil properties 

According to the logistic model in this study, the pres- 
ence and absence of earthworms can be a proper indica- 
tor of changes in soil physico-chemical characteristics. 
Soil salinity mostly results in the reduction of growth, 
community sizes and macro and microorganism activ- 
ity in soil (Owojori et al. 2009). This is due to the fact 
that soil salinity directly affects nutrition (Nemati et al. 
2011). Studies showed that earthworms can be an indi- 
cator of salinity stress in soils (Owojori and Reinecke 

2009). High contents of nutrient elements in Un and 
TCM sites and earthworm high abundance showed that 
earthworms can be a suitable indicator for soil fertility 
in Zagros forest ecosystem. Earthworms’ existence is 
related to accessible organic matter in soil (Poier and 
Richter 1992). Nurhidayati (2011) investigated the ef- 
fects of long-term forest management on earthworm 
diversity and density. The results showed that with 
25% increase of organic carbon, density and biomass of 
earthworms increased up to 64% and 83%, respectively. 
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Furthermore, with an increase of soil total N to 25% 
, density and biomass of earthworms increased up to 
79% and 75%, respectively. Curry (1998) stated that soil 
pH reduction is an effective factor in decline of earth- 
worm communities. Soil pH and soil physicochemi- 
cal properties affect the numbers and variety of soil 
organisms, and these factors are influenced by human 
activities such as land degradation (Brady and Weil 
2008). Soil aeration and oxygen supply are necessary 
for soil organisms to transform organic N to mineral 
N (NH4-N and N03-N) (Brady and Weil 2008). Earth- 
worms prefer habitats with proper aeration and their 
activities in soil lead to formation of horizontal and 
vertical tunnels, which improve soil aeration. In fact, 
earthworms are important as organisms which interact 
with microorganisms and influence primary produc- 
tion, decomposition and nitrogen cycling (Nadiri et al. 
2011). In the present study, both the logistic model and 
PCA proved that there was a good relationship between 
earthworm community and soil aeration, OC as well 
as mineral N (NH 4 -N and N0 3 -N). When compared 
to other sites studied, soil structure and BD in Un and 
TCM were more suitable, which can be related to high 
abundance of earthworms in these sites. The change of 
soil structure into granular shape is due to earthworms, 
who also crush plant debris, change its positioning and 
push forward organic matter parts enhancing formation 
of aggregates so that it can be considerably influences 
organic matter dynamics and soil chemical properties 
(Ladd et al. 1993). Based on multiple regression analy- 
sis, Tondoh et al. (2007) found that OC and soil pH were 
potential criteria for changes in earthworm abundance. 
Earthworms respond to soil degradation - equally to 
soil physical and chemical changes (Linden et al. 1994). 
Analyses of soil quality factors showed that long-term 
planting of tea and destruction of site caused a decrease 
of organic matter, N, P, K, porosity, soil acidity, and an 
increase of BD and soil compaction. Earthworm abun- 
dance as a biological index declined significantly due to 
long-term planting of tea (Van dang 2007). 

Along with the results of this research, it can be 
stated that human activities caused forest destruction 
and changes in physical and chemical properties of soil 
in the Zagros region. Abundance and presence and ab- 
sence of earthworms could be a suitable indicator for 
forest management. Hence, earthworms as such indictor 
can be used to monitor susceptible sites as Zagros, eas- 



ily and less time-consuming. Based on the results of 
the present study, it can be declared that unlike other 
soil characteristics, suchlike salinity and acidity (qual- 
ity indices of the forest site according to other studies), 
earthworms as a part of soil phone deliver more various 
physical and chemical properties simultaneously, and 
they might be a proper index in diagnosing soil quality 
in flam oak forest. 

Conclusion 

Management approaches should be evaluated with ap- 
propriate measures to determine how a given area re- 
sponds to management. In the present study, we evalu- 
ated the potential of earthworms as an indicator for dif- 
ferent forest management types and human disturbanc- 
es in Zagros oak ( Quercus persica Jaub. and Spach) 
forests. The study area had different soil attributes due 
to disturbance and short- and long-term conservation 
management. Even with conservation management car- 
ried out in subsequent 5 years, FCM site demonstrated 
no desirable conditions. In contrast, after 20 years of 
conservation management, TCM site conditions were 
analogous to those observed in Un site. According to 
our results, earthworm, biomass and abundance showed 
high sensitivity to BD and EC increase. Our study 
showed that earthworm abundance and biomass could 
be a good indicator to evaluate different management 
types in forest sites with different conditions. 
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